Survival Probability of a Local Excitation in a Non-Markovian Environment: Survival 

Collapse, Zeno and Anti-Zeno effects. 
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The decay dynamics of a local excitation interacting with a non-Markovian environment, modeled 
by a semi-infinite tight-binding chain, is exactly evaluated. We identify distinctive regimes for the 
dynamics. Sequentially: (i) early quadratic decay of the initial-state survival probability, up to a 
spreading time ts, (ii) exponential decay described by a self-consistent Fermi Golden Rule, and (iii) 
asymptotic behavior governed by quantum diffusion through the return processes and leading to an 
inverse power law decay. At this last cross-over time a survival collapse becomes possible. This 
could reduce the survival probability by several orders of magnitude. The cross-overs times ts and 
in allow to assess the range of applicability of theFermi Golden Rule and give the conditions for the 
observation of the Zeno and Anti-Zeno effect. 
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I. INTRODUCTION 

The decay of an unstable local state is usually de- 
scribed, within a Markovian approximation, by an ex- 
ponential decay with a rate given by the Fermi Golden 
Rule (FGR). However, this description contains approx- 
imations that leave aside some intrinsically quantum be- 
haviors Indeed, works on models for nuclei, com- 
posite particles 0] and excited atoms predict that 
the exponential decay does not hold for very short and 
very long times, and this exponential decay may shows 
superimposed beats. 

In Ref. Q we presented an exactly-solvable model de- 
scribing the evolution of a surface excitation in a semi- 
infinite chain. Physical realizations of one-dimensional 
systems are provided by electron transport in superlat- 
tices [5] , discrete diffraction in photonic crystals [y| , and 
spin excitations in a chain of nuclear spins under an XY 
interaction [3] or under double quantum interaction 
Here, we quantify and interpret the short and long time 
limits, ts and t^j, of the FGR. We identify three time 
regimes for the survival probability Poo(0- Initially the 
decay is quadratic and it holds up to a time ts- From ts 
to tn it is exponential, and finally, for long times, it fol- 
lows a power law. The time ts gives an upper bound to 
the time interval at which repetitive projection measure- 
ments could lead to a Quantum Zeno Effect (0, [H). On 
the other hand, at tfi, a dip in Poo(0 of several orders of 
magnitude may occur. This survival collapse is identified 
with a destructive interference between the pure survival 
amplitude, i.e., an exponential decay amplitude, and a 
return amplitude, which is usually neglected because it 
arises from memory effects in the environment. This de- 
structive interference can be used to obtain an anti-Zeno 
effect where the decay rate is strongly enhanced by 



repeated projective measurements with period tji. 



II. SURVIVAL PROBABILITY 

The evolution of a state |0) weakly coupled to a set of 
states which defines the "environment" , is described by 
the survival probability 



Poo (t) = mexp[-int/h]\O)0it)\' 



_de 
2ttH 



G^o(^)exp[-i£i/^] 



/OO 
de A^o (e) cxp[-iet/h] 
-OO 



(1) 

(2) 
(3) 

(4) 



'Electronic address: lrufeil@famaf.unc.edu.arl 
t Electronic address: lhoracio@famaf.unc.edu.arl 



where (t) is the retarded single particle Green's func- 
tion and A^o i^) is the Local Density of States (LDoS). 
This is evaluated expanding the initial condition in the 
eigenstates |fc) of the Hamiltonian 7i, or by using the 
energy representation of the Green's function Gqq (e), 

7Vo(e)^5]j(0|fc) 1^5 (5) 
--l/7rImG?o(e). (6) 

If the spectrum is bounded, Eq. ^ can be calculated 
using the residue theorem with the path shown in Fig. [1] 
Resonances appear as poles of the analytical continua- 
tion No{z) = No{e + ie') in the lower complex semi- plane. 
A well defined resonance appears when an initially un- 
perturbed state of energy Eq — (0|7i|0), far enough from 
the band-edge, is weakly coupled to a continuum, i.e., 
the expansion of |0) in terms of the eigenstates has a 
small breath Fq around an energy Sr = £o + Aq , where 
Ao = A{e — Er) is a small shift due to the interaction. 
This condition excludes out-of-band resonances, virtual 
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FIG. 1: (color online) Local Densities of States (LDoS) in 
the complex plane z — e + ie' . sl and eu are the lower 
and upper band-edges, respectively. The solid line is A'^o (e) 
for a semi-infinite chain with a surface impurity. The dotted 
hue is iVf ' (e) for a semi-infinite homogeneous chain, and in 
the inset is shown the lower band-edge of A^o (e). The pole 
appears in e^— iFo. The integration path is shown with dashed 
lines; it consist of four straight lines and two arcs (which avoid 
the band-edges singularities). 



states and localized eigenstates Then, 



SC-FGR 



(7) 



return correction from quantum diffusion 



where a = 27ri \im.z^er-iTa [{'^ — £r + iro) Nq{z)] is the 
pole residue. Pqo presents two separate contributions 
for the decay. The first term (the pole contribution) of 
Eq. ([7]) supersedes the usual FGR approximation since it 



has a pre-exponential factor [A 



> 1) and an exact 



rate of decay Fq, i.e., this result is a self- consistent Fermi 
Golden Rule (SC-FGR). This term is the dominant one 
for a wide range of times, leading to 



Poo (t) w Aexp(-2Fot/;i) 



(8) 



By analogy with the self-diffusion process in a classical 
Markov chain, the exponential in Eq. ([7]) is identified 
with a pure survival amplitude. Within the same anal- 
ogy, the second term (the integration path contribution) 
will be called return amplitude. The "quantum diffu- 
sion'^ described by this term dominates for long times 
and brings out the details of the spectral structure of the 
environment. 



A. Short time regime 

The second term of Eq. ([7]) is fundamental for the 
normalization at very short times. Both terms combine 
to provide the initial quadratic decay required by the 
perturbation theory 

Poo(^)-l-((£-£,.)')^^„^Vft' + ••■ , (9) 

where ((e — £r)^)^^ is the second moment of the LDoS 
Nq (e). This expansion holds up to the characteristic time 
ts- Let us consider a single state of energy £o coupled by 
Vb.j to an environment defined by N states of energy £j 
spread over a bandwidth B, as shows Fig.(I2]-a). 



(a) 



(b) 
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FIG. 2; (color online) (a) Single state of energy eo coupled 
by Voj to N states of energy Sj, spread over a bandwidth B 
with mean-level space of A « 2B/N, that defines the envi- 
ronment, (b) Equivalent semi-infinite chain with eo,Vo,ei, 
defined in text. 

A semi-infinite chain can be obtained from this system 
by using the recursion method [l^ , a variant of Lanczos 
tridiagonalization scheme, as is shown in Fig.([2]-b). The 
first two states are: 



eo = eo; |0) = |0) , 



(10) 



£1 = ^ - Sk^, (11) 



Vr 







Vo 



Vo 



(12) 



Here, the local second moment of the Hamiltonian is Vq , 
leading to 

Poo {t) = i- vftyh^ + ■■■ (13) 

There is a simple expression extrapolating Eqs. and 
®, El: 

Poo {t) « exp[(l - Vl + (t/tsr) ^^ots/h], (14) 

ts^h^. (15) 



with. 
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This yields Eq. ((13)) for t <C at the lowest order. In 
contrast, for t > is it yields the SC-FGR of Eq. with 
A « ex.p{2Tl/V^) valid for |K| < B. Therefore, as was 
remarked by Pascazio et al. [io|, the upper limit for the 
quadratic behavior is not H/Vq, as one might expect, but 
rather the much shorter time ts- A useful interpretation 
of ts can be drawn from the Green's function |14| : 



1 



e-eo 



R{0) 



(16) 



where G—"^ (e) corresponds to a semi-infinite chain in 
absence of Oth-site. Taking G^^°^ (e = £o) Eq. (fT6|l gives 
the FOR: Tfgr ~ nV^N^ (eo). Replacing it in Eq. 
(fT5|. we get 



ts 



hnNf (eo) 



(17) 



only determined by N-*^ (eo), the LDoS at the Ist-site of 
the unperturbed environment, evaluated at £o- In turns, 
SA^-"-* (eo) represents the time scale [l^ at which an ex- 
citation built from the decay, decays into the rest of the 
environment. Therefore, the return to the Oth-site, re- 
quired to build up the quadratic decay, becomes less ap- 
preciable than the escape towards the chain, leading to 
the fast exponential decay of the survival probability. 



B. Long time regime: Survival collapse 

For long times, only small values of e' contribute to the 
integral of the second term in Eq.([7]). This restricts the 
integration to a range near the band-edges. Then, taking 
into account Eq. ^ and performing the Fourier transform 
retaining only the Van Hove singularities at these edges, 
we get the power law decay at long times. This second 
term dominates Poo (t) because its decay is slower than 
the exponential one. The relative participation on the 
LDoS at each edge is /3 = [(e^-eL)^+rg]/[(et/-er)^+r§]. 
Collecting both edge contributions gives 

^^00(0 « [l + - 2(3cos{Bt/n)] 



— le 



(18) 



This means that the long time behavior is just the power 
law multiplied by a factor containing a modulation with 
frequency B/h. Eq. ([7]) shows that the survival am- 
plitude of the local excitation recognizes two alternative 
pathways: the pure survival (pole contribution), and the 
returning pathways where the excitation has decayed and 
explored the environment. Then, there is an interference 
term that becomes important when both amplitudes are 
of the same order. It is precisely at this cross-over time 
tfi between the exponential regime and the power law 



regime when the interference term can produce a sur- 
vival collapse, i.e., Pqc {t) nearly cancels out. This effect 
is seen as a pronounced dip in Fig. [31 

We also note that if the unperturbed energy state eo 
is exactly at the center of the band, [3—1, the pure 
return probability presents periodicals zeros barely com- 
pensated by the small pure survival probability. This 
should not be confused with the survival collapse dis- 
cussed above, which may yield an exact zero in Pqo (^)- 



III. SEMI-INFINITE CHAIN: EXACT 
SOLUTION 

Let us focus on a tight-binding Hamiltonian shown in 
Fig. (PJ-b) with hoppings Vq = Vo, Vjj+i = V and site 
energies eo and ej = 2V for j > 0: 

7^ = |0) eo (0|-(|0)Vb(l|+c.c.) 

+ (|n) 2V {n\ - \n) V {n + 1\ + c.c.) . (19) 

n 

This defines a continuous spectrum [el — 0,eu — 4F — 
B] and a well defined resonance for Vq ^ V. We first 
summarize the results in Ref. The LDoS factorizes 
as a pure Lorentzian around Sr ± iFg, and n['^\s): 



= 1^7 "r^-Ti^^^^)' (20) 

1 c (er - e) + 



with 



N[^\s)^'-^0i\e-2V\);Tie) = ^ 



e^ = eo-^Ao; Aq = '° 



Vq^ So - 2V 



(21) 
(22) 



r, 



r,; r, = dv^-Vo^ 



eo-2V 



(23) 



The solution of Eq. (U) results in: 



1 - {Vot/ny , t<ts 

Poo(t) w { Aexp{-2Tot/h), ts < t < tn 

C[l ^ j^sm{Bt/h)][h/{r{er)t)f,tn < t 

(24) 

with (3 as defined in section 2 and 



A = 



./elV^,^{B-erf+Tl 
V^^VV{erf{l+p^) 



(25) 
(26) 



Notice that the cubic power law decay at long times fol- 
lows from the -^e dependence of the LDoS near the band 
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edge (see inset of Fig. Fig. ^ shows Eq. for 

Vq/V = 0.4 and Sq/V = 1. One obtain an alternative 
representation of Eq.© by introdu cing an effective de- 
cay rate TcS (t) — —h/{2t)\nPoQ(t) [10'| whose deviation 
from Fq is a signature of non-exponential decay. This 
is shown in the inset of Fig. There, the survival 

collapse is a pronounced peak in Fcff (t). 




FIG. 3: (color online) Survival probability for a semi-infinite 
chain with eo/V = 1, Vo/V = 0.4, that leads to a resonance 
in Er/V = 0.9, Fo/V = 0.14. The inset shows Fe// (t) (solid 
line) and Fq (dashed line). 



and Fcff (t) for Eq/V = 1.8 and Vq/V = 0.77, which yields 
to tu « 6.8 [h/V]. In this case, recursive projective mea- 
surement at a time interval « ta can make the survival 
probability much smaller. Then the survival collapse en- 
ables an acceleration of the decay induced by repetitive 
observations, i.e., an anti-Zcno effect [l0|. 




FIG. 4: (color online) Survival probability for a semi-infinite 
chain with eo/V = 1.8, Vo/V = 0.77. The inset shows Feff (t) 
(solid line) and Fo (dashed line). 



In order to obtain the characteristic time in the 
weak coupling limit we solve iteratively the equality be- 
tween the exponential and the power law decay (aver- 
aged in a period), starting with h/2TQ. Since, for Sq 
close to the center of the band \J A/C ~ \/327ry/Fo and 
F {sr) ~ y ^ we obtain 

where oi , a2 > 1 are constants that depend on the Van 
Hove singularity (e) ~ (e — ei,Y ^-nd other details of 
the model. For a semi-infinite chain ai = + 2 = 5/2 
and ai = s/Ftt ^ 1.6. By choosing the parameters Vo,eo 
as above, this characteristic time results in k, 41 
\h/V\^ which is somewhat smaller than the exact time 
tji K, 62 ^/V\. Just the next order of iteration gives 
a much better approximation ^ 67 Also, by 

using these parameters, Eq. (fTT]) results in ts ~ 0.8 [h/V] 
which is a good bound for the short time scale. 

In the range of quadratic decay, recursive projective 
measurement of state |0) at a time interval would pro- 
duce a deceleration of the decay, i.e., a Quantum Zeno 
effect (see, for example. Fig. 1 in Our results above 

provides a convenient upper bound, < ts, for this time 
scale. 

The survival collapse can also occur at the strong 
coupling limit at the cross-over between the short time 
regime and the power law decay. Fig. (jH) shows Pqo (t) 



IV. CONCLUSIONS 

We studied the dynamics of a local excitation in a 
system in which full memory effects at the environment 
are included. We obtain the time limits where the non- 
exponential behavior of the survival probability shows 
up. The evolution starts with the expected quadratic de- 
cay, which holds up to a time ^5 (Eq. (fTT]) ) determined 
by the density of the first state of the environment in ab- 
sence of the initial state. This time gives an upper bound 
to the interval at which repetitive projection measure- 
ments leads to a Quantum Zeno effect. In the weak cou- 
pling limit the decay follows the usual FOR exponential, 
but with a corrected rate and a pre-exponential factor, 
i.e., the SC-FGR. At long times we get a power law decay 
controlled by non-Mar kovian return processes. We also 
visualized a survival collapse at time t^j (Eq. (P7)) ) as 
a destructive interference between the pure survival am- 
plitude and the return amplitude. This last arises from 
pathways that explore the environment before returning. 
Given that a survival collapse occurs, one can use repet- 
itive projective measurements with a period « t/j to 
achieve an anti-Zeno effect. 
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